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I. INTRODUCTION 



The improved characteristics of a modern missile make the interception process 
difficult: the reduced Radax Cross Section (RCS) makes it hard to detect by using 
a radar, while late emission leaves little time for intercepting and initiating any 
necessary defense and Electronic Counter Measures (ECM) activities. 

Conventional interception methods give only direction, and rough estimated 
range, and require a relatively long period of time. 

This work suggests a unique, non-conventional method for passive range and 
height estimation of a cruising missile, as well as any other microwave transmitter. 
The method is based on over sea multipath propagation, and uses five receiving 
antennas in a ladder configuration. For each receiving antenna there is a particu- 
lar multipath behavior, different than the other antennas. Comparing theoretical 
ratios of antenna pairs with actual received power ratios leads to an estimated 
location (range and height) of the transmitter. 

The motivation for the need of such a method is given in Chapter II. Chapter 
III describes the theory behind multipath propagation of electromagnetic waves, 
over different kinds of sea surface. It also explains in detail the model for rough 
sea which we used in this work. 

Chapter IV gives the basic theory of how the multipath propagation can be 
used for range estimation, and the conceptual estimation sequence. The detailed 
process, the assumptions being used during the research, and simulation results 
are given in Chapter V. 
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Additional parameters that were not taken into account in the process, as well 
as those only mentioned, are considered in Chapter VI, along with their influence 
on the experiment. The results are ainalyzed. 

Chapter VII explains the computer program, and shows the concept behind 
the software process. 

Chapter VIII summarizes and concludes the whole work, and gives some sug- 
gestions for future steps. 

The rough sea model which we used in this work is the “Ament Model.” Other 
familiar models are described in Appendix A. Listing of the computer program used 
is given in Appendix B. 
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II. NATURE OF THE PROBLEM 



A. GENERAL 

In the modern battlefield, where the electromagnetic complexity and density, 
as well as a missile’s speed are very high, the probability of intercept is low, so 
detecting and intercepting missions became hard to achieve. The policy is to start 
the transmission as close to the target as possible, zind to transmit for a minimum 
amount of time in order to be less “visible.” 

On the other hand, the radar cross section (RCS) of ships and especially 
missiles is being reduced, maJcing it difficult for a radar to detect them. 

Therefore, we see that intercepting a missile, especially an ultrasonic one, is a 
tough mission and might be achieved only during the Izist seconds, which might be 
too late. 

By using a conventional interceptor, even if it is a fast one (monopulse re- 
ceiver), only a direction can be found. Using the ship’s radar for range finding can 
be worthless, because of the small RCS of the missile, and at the same time will 
help the other side to detect the source. 

For these reasons, we prefer to find a missile (as well as other ships’) location 
(direction and range) by using high speed passive equipment. 

Two common methods for solving this problem are being used, and are listed 
below: 

1. Triangulation 

This method requires two independent interceptors on two different plat- 
forms sharing the information of the target direction. The location of the target 
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can be found on a map, as shown in Figure 2.1. Since each receiving antenna has 
its beamwidth, the position estimate is actually within a rhombus, and may not be 
as accurate as desired. The direction errors of the interceptors also add a particular 
level of inaccuracy to the estimation process. 

Other disadvantages of this method are the need for two ships to be in- 
volved, with all the logistics and synchronization, and obviously the time frame 
needed for this process, which maJce it a poor choice for missile defense. 




Figure 2.1: Triangulation Geometry. 

2. Auto-triangulation 

In order to estimate the location by using a single interceptor (single ship), 
we can make the triangulation by making sequential measurements, using the fact 
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that the interceptor itself is moving and chaxiging its place. Obviously, this method 
takes a long time for processing, and therefore is good only for stable targets (shore 
radar, for example). 

Actually, we cannot find the range of a moving target, especially a cruising 
missile, by using conventional methods, and so we need to have a special way for 
estimating the transmitter location. 

The method suggested in this work uses the multipath propagation of a 
transmitted electromagnetic signal in order to estimate the transmitter range from 
the receiver, while the direction is still being found using any common method such 
as rotating antenna. Instantaneous Direction Finder, etc. 

The theory of a general multipath propagation over smooth and rough sea 
is described in the next chapter. 
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III. THEORY OF MULTIPATH PROPAGATION 



OF ELECTROMAGNETIC WAVES 



In a usual radar case, the multipath propagation generally causes tracking 
difficulties, especially over a smooth sea, which causes specular reflections. 

In this chapter we will review the multipath propagation theory, divided into 
smooth and rough terrains. 

A. GENERAL 

A transmitted signal from a radar antenna spreads as it is leaving the antenna 
surface. A direct path exists between the transmitter (radar) antenna and the 
target (receiver antenna in our case), as well as an indirect path along which the 
signal transmitted from the transmitter antenna hits the sea surface and is reflected 
to the receiver antenna. 

The signals arriving from these two paths create interference in the receiver, 
and the relative phase difference dictates if it is a creative or destructive inter- 
ference. The path difference depends on the direct range R, transmitter antenna 
height ht, and receiver emtenna height hr. Other factors have some influence as 
well, and will be discussed later in the next sections. Figure 3.1 shows the geometry 
of the multipath propagation over a smooth sea, assuming a flat earth. 

Since the sea state has significant contribution to the interference, a basic 
question that can be asked is how to distinguish between a smooth and rough 
surface (sea in our case). The answer for this question is not a straight-forward 
one, and, as a matter of fact, there is neither a single answer nor a defined border 



6 




IMAGE 
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Figure 3.1; Simplified Multipath Geometry Over Smooth Terrain and 
Flat Earth Assumption [From Ref. 1]. 

between the two cases. A method for having a quantitative solution for this problem 
was suggested by Lord Rayleigh, and is described herein [after Ref. 2]. 

Before starting the detailed analysis, we will give some definitions which are 
relevant to our discussion, referring to Figure 3.2: 




Figure 3.2: Geometry of Irregular Surface and Path Difference [From 
Ref. 3]. 
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• Wave Height A/i: The vertical distaxice between the wave’s trough and its 
crest. 

• Wavelength 1: The horizontal distajice between two adjacent crests (measured 
in the direction of wave propagation). 

• Wave Period T: Time interval between two adjacent crests passing the same 
fixed point relative to the groimd. 

Consider an irregular surface with irregtilarities of height Ah. If two parallel 
rays of the same wavelength A hit the surface with a grazing angle xf), as illustrated 
in Figure 3.2, the path difference between them will be Ar, where: 



Ar = 2A/isinV’. (3-1) 



Translating this into phase difference, we obtain: 



A<l> = kAr = ^Ar, 



and substituting (3.1) we have: 



47tA/i . , 

A<p = — - — smxp. 



(3.2) 



The phase difference A(f> can vary between 0 and tt, where a totally smooth 
surface will occur when A(f> = 0. Zero phase can be achieved when there are no 
irregularities, i.e.. Ah = 0, or when the grazing angle is zero. From the above, we 
can assume a smooth surface when: 



A 



0 or xj; 0. 



(3.3) 



When the phase difference between the direct and reflected signal components 
is equal to tt, there will be a destructive interference, and no received energy at all 
if there is a perfect reflection (|/9| = 1). 
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In order to choose a point where the surface becomes rough, we can arbitrarily 
choose any value between A(f> = 0 and A<f> = tt. The most popular value is half the 
way (A<j> = 7t/2), as was chosen by Rayleigh, and hence, the “Rayleigh criterion” 
for a smooth surface is: 



Ah < 



(3.4) 



8sin^ 

Other values have been chosen as well [Ref. 2,3,4,5,6]. 

From observing Equation (3.2), we can learn that the phase behavior is a 
periodic one, hence the interference will have maxima and minima in sequence. 

A typical received signal in the presence of multipath propagation over a 
smooth sea is shown in Figure 3.3. The amplitude decay is due to the spheri- 
cal spreading of the received signal (1/R^). 



MUl:iPA1H effect - SEA STATE 0 
RX. ANT. HEIGHT: IfiM 




Figure 3.3: Multipath Propagation. 
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B. SMOOTH SEA 



Looking at the multipath propagation more carefully, we should notice that 
there are additional phase and amplitude changes in the reflected wave due to the 
reflection itself, as well as amplitude difference between the reflected wave and the 
direct one, since the propagation path lengths are not the same. 

The amplitude difference due to the difference in path lengths is very small, 
and hence can be neglected, since the paths differ only by a short distance compared 
to the total path length. [Ref. 7] 

The phase and the amplitude changes due to the reflection itself depend on 
the grazing angle, frequency, and polarization. 

Figures 3.4 and 3.5 illustrate the magnitude and phase (respectively) for verti- 
cal polarization, while Figure 3.6 shows the magnitude of the reflection coefficient 
for horizontal polarization. 

Note the small changes in the value of the reflection coefficient for horizontal 
polarization, which means that amplitude changes due to reflection can be neglected 
for small grazing angles (assuming reflection coefficient magnitude = 1 ). The phase 
difference in horizontal polarization can be considered as 180°, since for grazing 
angles of 0° to 90° it changes only up to 4°. [Ref. 7] 

Assuming a flat earth and a smooth terrain (reflection coefficient = — 1), the 
ratio between the received power in the receiving antenna in the presence of mul- 
tipath, compared to a free space (for small grazing angles) is [Ref. 3,9], where; 

=4sm ^ (t^ 0), (3.5) 

using the same notation as in Figure 3.1. The shape of rp' is shown in Figure 3.3 
for a typical case. 
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Figure 3.4: Reflection Coefficient in Vertical Polarization 
[From Ref. 7]. 

In order to be more precise, we need to taJce into account the fact that the 
earth is not flat. This fact will cause some changes in the geometry, and hence will 
influence the signal amplitude. 

Actually, the reflection coefficient depends on three main parameters: 



1) Plane smooth surface reflection coefficient p. 

2) Earth curvature, Z>,which reduces the reflection (in some cases, depending on 
the geometry). 

3) Surface roughness Rn. 

The actual reflection coefficient, hence, is [Ref. 6]: 



p = p - D ■ Rn. 



(3.6) 
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Figure 3.5; Phase Difference in Vertical Polarization [From Ref. 7]. 

The value of R„ is between 0 and 1, depending on the terrain, and it will be 
discussed in detail in the next section. 

The earth curvature factor D (Divergence factor) has also values between 0 
and 1, where in the flat earth case it will be unity. 

Practically, in the scope of this work, we can use the flat earth assumption, 
since the geometry of the problem (transmitter height less than about 100 meter, 
receiver antenna height at about the same and both heights are much smaller than 
the distance in between) leads to small grazing angle and allows the assumption of 
parallel direct and reflected rays [Ref. 7]. 

Even though the transmitter is considered a point source, the reflections from 
the sea come from a particular zone, rather than in a single point. We shall look 
for all the points (loci), from which the reflected signal will reach the receiving 
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Figure 3.6; Reflection Coefficient in Horizontal Polarization 
[From Ref. 8]. 




Figure 3.7: First Fresnel Zone [From Ref. 2]. 

antenna with the same phase shift, relative to the direct radiation, as illustrated 
in Figure 3.7. The loci obey Equation (3.7). 
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R\ + i?2 Rd — 



(3.7) 



or after arranging: 

Ri "f* i ?2 ~ R-d “I" = Constant. (3*^) 

\ 

Equation (3.8) describes ellipsoids, so that the reflection zones are actually 
ellipsoids which form the Fresnel zones (Fresnel ellipsoids). Near the wavelength 
that has the most interest to us, 3 cm, (or frequencies around 10 Ghz), the area 
of the first Fresnel zone is about 2000 square meters [Ref. 2,3,6]. It is common to 
take the whole flrst Fresnel zone as the reflection area. 

C. ROUGH SEA 

The smooth sea case is generally rare, while some kind of roughness is more 
common, producing complications to the problem by adding more parameters, 
some of them arbitrary. Due to these facts, the rough sea case is difficult to 
analyze and predict. 

A more accurate solution leads to more complications and difficulties, so in 
practice we try to simplify and make assumptions in order to get a solvable case. 
Still we need to be very careful so that the result will approximate reality even 
after using these assumptions. 

In general, we can say that the first problem is to find a correct compromise 
between excess rigor which may yield no practical result, and too many simplifica- 
tions which may lead to incorrect results and hence wrong conclusions [Ref. 3). 

In practice, the common general assumptions which take place in the rough 
sea case are: 

• Scattering elements are much larger than the propagated wavelength. 
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• Only the far field is calculated, 

• Diffraction and reflections between waves axe neglected, 

• A particular model for the roughness shape is assumed (depending on the 
case), and average wave height is assumed, and, 

• Constant sea state (roughness) is assumed for the whole area of interest. 

1. Reflection types 

At this point we will distinguish between two kinds of reflections: 

a. Specular reflection 

Similar to the smooth surface (as in a mirror case), specular reflec- 
tion can be characterized as obeying the classical optics rules, as well as having 
a coherent phase. The direct source for this kind of reflection is the first Fresnel 
zone (Fresnel ellipse) radiated points. The phases of the different points (from the 
receiver point of view) are approximately equal, hence the reflection is coherent, 
as mentioned previously. The specular reflection fluctuations are relatively small 
in amplitude. 

b. DiS'use reflection (scattering) 

Caused due to irregularities in the surface, scattering cannot be pre- 
dicted from the phase and directivity points of view. This scattering will occur 
over larger areas than the specular reflection (first Fresnel zone) and is incoherent 
in phase. The diffuse scattering has large amplitude fluctuations. We can also say 
that the smooth sea analysis is deterministic, while the rough sea analysis (both 
specular and diffuse reflections) is stochastic in nature [Ref. 5]. 

In general, there will be both reflection components, since a surface is 
neither totally smooth, nor extremely rough. The two components are illustrated 
in Figure 3.8 [Ref. 2,7,10]. 
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The specular and diffuse reflection coefficient components will be de- 
noted as p 3 and respectively. The effective reflection coefficient will be a mul- 
tiplication of the smooth terrain coefficient po» the relevant specular or diffuse 
coefficient. Hence, for the specular component, pj, the reflection coefficient will be: 

p = Po • P*, (3.9) 

and obviously will be smaller than that of a smooth sea po (since each coefficient 
is less than unity). 

The diffuse reflection will have a coefficient of 

p = po-pd, (3.10) 



again, smaller than po- 




Figure 3.8: Illustration of Diffuse and Specular Reflections 
[From Ref. 10]. 
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The diffuse reflection coefficient, due to its random nature, is found 
to have a Rayleigh distribution probability function for the amplitude, as given in 
Equation (3.11), while the phase has a uniform distribution between 0 and 27 t. 
The diffuse reflection coefficient amplitude has the form: 

Pd = v/2. r • exp(jVa), (3.11) 



where: 



cTa is the RMS wave slope (unitless), as shown in Figure 3.2, 
r = 0.77 • [1 — exp(— 47 t^)] • exp(— 4.73 • 5), and, (3.12) 

S = (3.13) 

A 

where ah is the standard deviation of the wave height, related to the peak to trough 
height. [Ref. 5] 

Assuming a smooth terrain, there will be a center point of reflection. 
The diffused reflection comes from the whole area surrounding this point, the glis- 
tening surface [Ref. 2,11], which depends on the observation angle from the trans- 
mitter to the receiver {Et in Figure 3.1), and on the RMS wave slope (<7a). 

A typical RMS slope, a a, is about 0.2 radians, and it depends on the 
wind velocity. The amplitude of the reflection coefficient has an extreme of 0.4, 
which rarely occurs, while a typical value is 0.35. [Ref. 11] 

We should notice that for a flat earth assumption, the grazing angle 
xj} is the same as the observation angle Et from the transmitter to the receiver 
(see Figure 3.1). The idea of the different reflection areas due to observation angle 
Ef, and sea state (roughness) is demonstrated in Figure 3.9. In this Figure, a a 
represents the rough surface RMS slope. 
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The specular reflection coefficient has been analyzed by many, and 
there are quite a few models for this component. Each model uses different as- 
sumptions for the computations, and therefore there are (in some cases) different 
results. The conclusion from this is that the results are highly dependent on the 
roughness type which is assumed. 
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Figure 3.9: Reflection Areas in Different Roughness and Observation 
Angle, Et> (a) Specular reflection, ~ l,p = po» (b) Slight 
roughness, (cr/i/A) • sin Ft ~ 0.06, (c) Rough Surface [From 
Ref. 11]. 



In order to be as close to the experimental data with a generic formula- 
tion as possible, we can, in general, define the reflection coefficient, p, as [Ref. 12]: 

OO 

P— p{4>) cos <f> d(f> , (3.14) 

— OO 
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where: 



^ is the phase of each particular element of the reflected wave in relation to a 
median plane surface, and, 

p{(j>) is the probability distribution of all the phases (f> along the path. 

2. Profile Types 

Using Equation (3.14), we find that the common profile models for rough 
terrain are [Ref. 2,12]: 

a. Sawtooth profile 

The phase (j) varies uniformly between maximum values of 
the probability distribution therefore is: 

= 51 -' 

while the reflection coefficient is: 

p = = sinc<?i,„. (3.16) 

YTTl 

The idea here is that the wave height is uniform, but wave slopes are arbitrary. 

b. Sinewave profile 

The phase <f> behaves like a sine wave with maximum amplitude of 
In this case, the probability distribution is: 

1 

2lZ(j)m COS U ’ 

where u is an arbitrary parameter, and: 

P ~ <^0 (^^m )? 

where Jq is the Bessel function of zero order. 





19 



c. Random profile 



The phase (f> has normal (gaussian) probability distribution with stan- 
dard deviation of cr. Here, 

exp ) 






\/^ ■ 



(3.18) 



and the reflection coefficient in this case is 



p = exp-(crV2)- (3.19) 

The latter profile is similar, in a sense, to the Ament model, which will be discussed 
later in this chapter. The three profiles are drawn in Figure 3.10 [Ref. 12]. 




IT iir 3T 



Figure 3.10: Reflection Coefficients of Common Models for Irregular 
Terrain [From Ref. 12]. 
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A further discussion on other models is given in Appendix A. 

The most general and the most popular model is the one by W. S. 
Ament.* This model assumes a normal (gaussian) distribution, and it is popular 
since it does not include any discontinuities as other distributions do, but still 
meets the actual (experimental) data. For these reasons it seems to be the most 
realistic model [Ref. 2]. 

The Ament model is discussed and analyzed in many books and arti- 
cles, and from now on we will concentrate especially on it [Ref. 2, 3, 5, 6, 7, 9, 11, 
13, 14, 15]. 

D. AMENT MODEL FOR IRREGULAR TERRAIN 



The Ament model assumes that the wave height is normally (gaussian) dis- 
tributed, and it does not take into account the shadowing effect, where parts of the 
surface which are hidden by higher parts are not being illuminated by the incident 
wave. The shadowing effect takes place at very small grazing angles, while for 
higher angles it has no effect. The geometry is shown in Figure 3.11. [Ref. 13] 
Figure 3.12 shows the influence of the shadowing effect on the wave height 
standard deviation, related to the grazing angle [Ref. 13]. 

The specular component of the reflection coefficient in this model is: 



Pa = exp -2 



( 2-K(Jh sinxp' 
A ^ 



(3.20) 



where: 

ah is the RMS deviation of the waves height, 
tp is the grazing angle (refer to Figure 3.1), and, 

* W. S. Ament, “Toward a Theory of Reflection by a Rough Surface,” 
Proc. I.R.E., vol. 41, p. 142, January 1953. 
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Figure 3.11: Shadowing Effect Geometry [From Ref. 13]. 




Figure 3.12: Shadowing Effect on Wave Height Standard Deviation 
[From Ref. 13]. 

A is the incident wavelength. 

Figure 3.13 shows the specular reflection coefficient component according to 
Ament model [Ref. 7]. 

The parameter ah is proportional to the sea state, and is given with other 
relevant parameters in Table 3.1 [after Ref. 2,11]. 
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Figure 3.13: Rough Surface Reflection Coefficient According to Anient 
Model. 



TABLE 3.1: WAVE HEIGHT PARAMETERS AND CRITICAL 
ANGLES FOR SPECULAR REFLECTION, VERSUS SEA STATE. 



sea state 
number 


description 
of sea 


wave 

height (m) 


RMS height 
(Th (m) 


critical an; 
3 GHz 


?le (°), lAmax 
10 GHz 


0 


calm 


0 


0 






1 


smooth 


0-0.3 


0-0.065 


> 6 


> 1.8 


2 


slight 


0.3-1 


0.065-0.21 


1.8-6 


0.5-1.8 


3 


moderate 


1-1.5 


0.21-0.32 


1.2-1.8 


0.3-0.5 


4 


rough 


1.5-2.5 


0.32-0.54 


0.7-1.2 


0.2-0.3 


5 


very rough 


2.5-4 


0.54-0.86 


0.4-0.7 


0.12-0.2 


6 


high 


4-6 


0:86-1.3 


0.3-0.4 


0.09-0.12 


7 


very high 


> 6 


>1.3 


< 0.3 


< 0.09 
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Figure 3.14: Scattering Regions Versus Roughness [From Ref. 11]. 



Using the Ament model, we can distinguish more carefully between rough and 
smooth sea, and check the conditions for diffuse and specular reflection. Recall 
(3.20) for convenience: 



Ps 



exp —2 



( 2'ir<Jh sint/» 
A . 



We can consider the terrain as a rough one if p, < 0.5, which is equivalent to: 



^ sin V’max ^ 0.0665. 
A 



(3.21) 



For a particular wavelength at a given sea state, there is a critical grazing 
angle V’max i calculated from (3.21), above for which the sea surface is taken as 
rough. This angle is shown in Table 3.1. 

The specular reflection coefficient component is significant until the grazing 
angle exceeds twice the critical angle. Beyond that point, the diffuse component 
becomes the significant one. The different regions of roughness type are illustrated 
in Figure 3.14. [Ref. 2,11,14] 
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Another way of looking at the specular (coherent) and diffuse (incoherent) 
scattering regions is given in Figure 3.15. Here we can see the boundaries between 
the scattering regions for different grazing angles. The graph in Figure 3.15 is given 
for 2 GHz. For higher frequencies, the graph will become more compressed, by the 
reciprocal of the frequency. [Ref. 15] 




Figure 3.15; Scattering Regions Envelopes at 2 GHz [From Ref. 15]. 

A similar graph for a constant wave height of 0.92 meter and different frequen- 
cies is given in Figure 3.16. Here again, for higher wave heights, the envelopes will 
be compressed by the same factor by which the height is increased. [Ref. 15] 
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Figure 3,16: Scattering Regions Envelopes at Peak to Through Wave 
Height of 0.92 meter [From Ref. 15]. 
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IV. THEORY OF RANGE ESTIMATION 



USING MULTIPATH EFFECT 



A. THEORETICAL SOLUTION 

The multipath effect attenuation, as a function of range over a smooth sea, is 
well known, if both the transmitter and receiver height, as well as the frequency of 
the signal are given. This was described in the previous chapter. 

For a rough sea, the overall shape remains similar, except that the curve has 
more fluctuations. Still, by using a particular model, such as the Ament model in 
our case, we can deal with the multipath propagation as a known factor. 

Figures 4.1 and 4.2 show typical multipath propagation over a calm sea (state 
0) and over a very high sea (state 7), respectively. We can notice that in the rough 
sea case, the attenuation coefficient decreases arbitrarily in the short ranges, and 
has an arbitrary ripple over the whole range. 

From the above, the multipath attenuation can be written as a function of 
the range between the transmitter and the receiver R and of the transmitter height 
/it, for a given receiver height hr'. 






/ir=const. 



(4.1) 



The received signal is attenuated relative to free space. It is a combination of 
both the direct and indirect components. Hence, the received signal can be written 
as: 

Sr = Sd-rj\ (4.2) 
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where: 

Sr is the received signal power, 

Sd is the direct component of the signal power, and, 

77^ is the multipath power attenuation factor, given by Equation ( 3 . 5 ). 

Since, as mentioned before, we cannot specify the direct component Sd, we 
need to eliminate it in order to solve the problem. A simple way for doing that is 
by using a pair of receiving antennas, one above the other, so that the range R will 
remain the szime, while their heights have no practical influence on the direct signal 
Eimplitude. We should also remember that the transmitter height, ht, is constant 
for both receiving antennas. 

If we use the ratio of the output powers from the two antennas, 

(■^r)i _ Sd'Til _tjI . . 

( 5.)2 Sd^vl vV 

we will be dependent only on the geometry of the receiver and the transmitter, 
without caring about the direct signals’ amplitude. 

Referring to Equation ( 4 . 1 ), we notice that we still have two unknowns, R and 
ht. Using the fact that we have two receiving antennas, while the range R and the 
transmitter height ht axe the same, we need to And the correct ratio 771/772 which 
will give the right combination of range and height {R and ht). Figure 4.3 shows 
the geometry of the problem. A more detailed explanation of the sequence is given 
below. 

Let us assume, for simplicity, that the multipath propagation curve is a linear 
function (we shall return to the actual case later on), hence it behaves like an 
isomorphic function. 
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Figure 4.1: Multipath Propagation Over a Calm Sea. 
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Figure 4.2: Multipath Propagation Over a Very High Sea. 



30 



(KM) 



TX. 




Figure 4.3: Range Estimation Basic Geometry. 

The power ratio of the outputs of the two receiving antennas can be written 
as: 

(5,)i S(R,h,)\K. Jhri\ , , 

(S.)2 nl /W*.)U,2 ^ ' 

SO it actually depends on the receiving antenna heights, which are known. 

Since we assumed that the problem has an isomorphic characteristic, a given 

ratio has only a single combination of range and transmitter height which fits this 

ratio. 

Keeping in mind those conclusions, the range estimation process will require: 

1) receiving a signal and analyzing its frequency and direction. 

2) using Equation (3.5), creating two (offline) tables (matrices) at that particular 
frequency, for which the coordinates are range and transmitter height, while 
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the contents axe the attenuation coefficients for each R and combination, 
while the parameter is the receiver height (hence there are two tables). 

3) calculating an (off line) look-up table for the same frequency, ranges and trans- 
mitter heights, which has the values of the ratio between the two tables created 
in step 2. 

4) calculating the ratio between the two received signals, which is similiar (see 
Equation 4.3) to the ratio of the attenuation coefficients. 

Now, all we have to do is to search in the look up table from step 3, and find 
the range and the transmitter height, so the last step will be: 

5) searching for the right combination of range and height. 

B. PRACTICAL SOLUTION (OVER SMOOTH SEA) 

Practically, the multipath propagation leads to a multi-valued curve, hence 
the problem is not isomorphic any more, but has an ambiguous solution. Since 
the function is not linear einy more, and involved with two unknown variables, 
we cannot use any analytical solution (like we could have used in the theoretical 
case), but only a numerical solution can take place. The idea is similar to the one 
described in the theoretical case with a pair of antennas. 

We saw before that two receiving antennas give a single ratio, hence we cannot 
solve the ambiguity problem (because of the periodic behavior of the multipath 
effect). 

Adding a third antenna will give a second relation, so we can now use the 
process described previously for the two pairs, with a few candidate combinations 
{R and ht), and by comparing their locations between the two lists of candidates, 
choosing that one which appears in both of the lists; all the other locations are 
false ones. 
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If the solution still is not a single one, we can add another antenna(s), until 
the comparison between all of the ratios leads to a single location. An important 
condition is that all of the receiving antennas be installed one above the other. 

A trade-off between the number of antennas and accuracy should be done: 
on one hand we want minimum ambiguous points (different points with the same 
ratios group) which dictates many antennas, but on the other hand we should use 
a realistic number of antennas. 

The influence of the different parameters on the multipath behavior can be 
seen in Figures 4.4 to 4.6. 

From observing these figures, we can conclude that: 

1) as the frequency increases, the curves become more complex (signal strength 
varies more rapidly with range), 

2) complexity increases with transmitter height, and, 

3) for increased receiving antenna height, the curves become more complex. This 
fact suggests keeping the antennas at lower height, but on the other hand, we 
should separate them enough to reduce the ambiguity. The receiving antennas 
in Figure 4.6 are 0.2 meters apart, so the curves are very similar to each other. 
In Figure 4.7, the antennas axe separated by 3 meters, and the curves are well 
distinguished. 



Note that for short range we have very high complexity, which will be treated later 



on. 



C. PRACTICAL SOLUTION (OVER ROUGH SEA) 

In the case of rough sea, the main idea remains the same, although the signal 
exhibits fading as a function of the sea state. As the receiving antenna is raised, 
the influence of the roughness becomes more important. Figure 4.8 shows the 
multipath propagation curves over sea state 7, for different antenna heights. It 
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is obvioiis from these graphs that the higher the antenna is, the more rapid the 
fading. 

The look-up tables should take in account also the sea roughness parameter, 
where in our case we choose the Ament model, as described in Chapter III. 

After considering the facts mentioned above, the estimation process remains 
the same as described in the previous sections. 
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Figure 4.4: Variation of Signal Strength with Range, for Different 
Transmitter Heights. 
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Figure 4.5: Variation of Signal Strength with Range, for Different 
Frequencies. 
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Figure 4.6: Variation of Signal Strength with Range, for Different 
Receiver Heights (Small Separation). 
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Figure 4.7: Variation of Signal Strength with Range, for Different 
Receiver Heights (High Separation). 
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Figure 4.8: Variation of Signal Strength with Range, for Very High Sea 
and Different Receiver Antenna Heights. 
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V. EXPERIMENTAL PROCEDURE AND RESULTS 



Now that we have presented the theory of how to use the multipath effect 
for estimating range, we can discuss how we applied that for solving the problem 
in our research. The assumptions, practical difficulties, and then some simulation 
results will be presented. 



A. ASSUMPTIONS 

Since the actual problem is a very complicated one, due to many surface 
parameters, where some of them are unknowns and/or very difficult to predict, we 
made some logical assumptions in order to obtain a solution. 

• Polarization: the seeker was assumed to have horizontal polarization, hence 
the reflection coefficient wzis taJcen as unity in amplitude, with a 180® phase 
shift. The differences for vertical polarization will be discussed in Chapter VI. 

• Sea state: the whole area of interest, i.e., the whole range between the trans- 
mitter and the receiver was assumed to have the same sea state condition. 

• Earth curvature: a plane earth was assumed. A spherical earth case is dis- 
cussed in Chapter VI. 

• Atmosphere: a homogeneous atmosphere was assumed. 

• Movement: the suggested method considers relative movement between the 
transmitter and the receiver. For homing missile, this is obvious. In other 
cases, such as estimating range between own ship and a stable transmitter 
(shore based radar, for example), we shall assume that the ship is moving 
relative to the radar. If there is no relative movement at all, the result can be 
incorrect. 

• Geometry: It was assumed that the antenna separation, from the point of 
reflection, was much larger than the signal wavelength [Ref. 2]. Since we 
are dealing with microwave signals (from 2 GHz at least), this assumption is 
valid. (If this condition is not true, the reflection coefficients have some more 
complexity.) 
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• Range interval: almost any range interval can be assumed, but note that 
for the long ranges we should check the geometry and see if the plane earth 
assumption is still applicable. Since we concentrated in this work on homing 
missiles, the seeker transmitting range is the one of importance, and hence a 
maximum range of 20 Km (about 11 NM) is considered, while a minimum of 
1 Km (about half a NM) assumed. 

B. ESTIMATION PROCESS 

With the assumptions of the previous section in mind, let us repeat briefly the 
range estimrtion process. We shall look more closely at the details shortly. 

1. General Discussion : 

A general solution requires the following steps: 

a. Frequency Measurement: 

Since the whole process is frequency dependent, this will be the very 
first step made by any EW receiver. 

b. Oflf Line Tables: 

In this step we build the theoretical tables for all ranges and heights of 
interest (i.e., the possible locations of the transmitter), for each receiving antenna 
height. The inputs for creating those tables are frequency, receiving antenna height, 
and sea state, which dictates the attenuation factor to be placed in the table. 

c. Attenuation Coefficient Ratio Tables: 

After we have the theoretical receiving values (from the previous step), 
we calculate (for each range and height combination) the ratios of the output signals 
for each receiving antenna pair. 

d. Amplitude Measurement; 

We must measure the received signal amplitude of each receiving an- 
tenna (we will discuss later the number of receiving antennas), and make a table of 
amplitude ratios. For example: suppose we have three receiving antennas, where 
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each one receives a signal with amplitude 5i,S2i -33, respectively. The ratio 
table will have two elements: 32/^1 and 53/31. 

e. Comparison: 

Now, that we have the measured signal ratio (from step d), and the 
theoretical one (step b), we need to compare - for each receiving antenna pair, the 
actual value with the olf line (theoretical) value, and record range-height combina- 
tions for which they axe equal. These range-height combinations represent possible 
transmitter locations. 

Since, as we mentioned in Chapter IV, the problem has an ambiguous 
solution, there will be several possible locations, even though only one of them is 
the correct one. The next step is to determine the correct location. 

f. Validation: 

In order to choose the correct location from all the candidate locations 
determined in step e, we will compare all the locations determined for the first pair 
of receiving antennas, and the second pair of receiving antennas. There is, ideally, 
only a single location, which should appear common to all the pairs. 

Practically, because of the high ambiguity level (especially for the rel- 
atively short ranges), and also because of approximations in the implementation, 
there will still be more than a single location. In order to resolve this problem, 
we shall use more receiving antennas, which provides more ratios. This allows the 
ambiguity to be reduced to an acceptable minimum. 

g. Output: 

The matrix location is translated into an estimate of transmitter range 
and height. This is accomplished by using the coordinates of the look up table. 



42 



2. Detailed Discussion 



The following section explains in detail every step listed before. 

a. Frequency Measurement: 

Since this requires an external system, we will not discuss it. Any EW 
receiver on board can be used for frequency measurement and direction finding. 
The frequency data is provided as an input to the range estimation process. 

b. Oflf Line Tables: 

This step uses the theoretical equations for calculating the multipath 
attenuation factor rj^ . For each receiving antenna we have a particular table, where 
the coordinates are transmitter height, ht (vertical axis), and range, R (horizontal 
axis). 

The height axis is divided into 36 “slices” with a 1 meter resolution, 
starting at a minimum height of 5 meters, and ending at a maximum of 40 meters. 
(These numbers are parameters, which can be changed easily.) 

The range axis runs from 1 Km to 20.5 Km, and was originally divided 
into 130 equal steps of 50 meters each. Since the ambiguity for short ranges is 
very high, this partition caused some problems in finding the correct location for 
short ranges. This suggested the need for higher resolution at short range. The 
chosen partition is 130 steps of 50 meters each, up to the range of 7.5 Km, and 
then reduced resolution of 125 steps, 100 meter each, until the maximum range. 
(Again, all these numbers are parameters which can be changed.) 

The multipath attenuation factor is calculated by using the following 

formulae: 



The basic attenuation, for a smooth sea case is [Ref. 3,9]: 



r) =4 sin 



. 2 2TTht{ht + hj.) 



\R 



(V’ ^ 0), 



(5.1) 
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where: ht is the transmitting antenna height, 
hr is the receiver antenna height, 

A is the signal wavelength, 

R is the range between transmitter and receiver, and, 
t/> is the grazing angle. 



coefficient, p, which is different in various models. In our case we used, as described 
in Chapter III, the Ament model: 



where: ah is the RMS deviation of the waves’ height, 
tp is the grazing angle, and, 

A is the signal wavelength. 



Therefore, in a rough sea case, the multipath attenuation factor will 



In the case of a rough sea, we have to multiply by the reflection 




(5.2) 



be: 




(5.3) 



The grazing angle tp in equation (5.2) is a derivative of the geometry. 



It depends on transmitter height, receiver height and range. Figure 5.1 shows the 
geometry from which we calculate rp. 



We shall now derive the grazing angle tp [According to Ref. 7]. From 
Figure 5.1 we see that 




We define the range difference between the indirect path and the direct 



path as <5, where: 



6 — R\ “I" /?2 — Rd- 



(5.5) 
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Figure 5.1; Flat Earth Multipath Geometry [From Ref. 7]. 

Looking at the image of the transmit antenna, we see that 



i?i + R 2 — \/ + B?' . (^-6) 

Substituting Equations (5.4) and (5.6) into Equation (5.5), we obtain: 
8 = ^{ht + KY + R2 _ ^i^ht - hrf + R } . (5.7) 



Let us assume that the range R is much greater than the antenna 
heights, and also R} (h< + hrY . In this case, we can approximate Equations 
(5.4) and (5.6) as: 



- KY + R} 



R 



{ht + hrY 
^ 2R? 



(5.8) 
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and, 



Rl + i?2 — {hr + ht')^ ~ i? 



1 + 



{ht - 

2B? 



(5.9) 



respectively. 

Substituting Equations (5.8) and (5.9) into Equation (5.5), we obtain, 
after arranging: 



{ht^-hrY -{ht-hrY 2hthr 



2R R 

Observing again the imaged antenna, the angle $r is given by 



(5.10) 



Or = sin 



— 1 hi + hr 



(5.11) 



and we notice that it equals (interchangeable angles) the grazing angle, t/>. Using 
Equation (5.5), we finally find the grazing angle to be: 



, . -I hi + hr 

xp = sin r. 

Rd ^ 



(5.12) 



Since we know all the parameters from solving previous equations, we 
can now find xp. (Actually, we do not know the range R, but we should remember 
that in this step we are building tables, so the range is being used as a parameter.) 

The only variable still missing for solving Equation (5.2), (the reflection 
coefficient), is the wave height RMS deviation cr^. This value depends on the sea 
state, and is given in tables in the literature. Table 3.1 shows this parameter (as 
well as some others), as a function of sea state. 

The propagation attenuation factor ( in the one way case) also 

has its impact on the received signal. Since we deal with ratios, this factor cancels 
itself, and therefore had not been taken into account during the calculations. 
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c. Attenuation Coefficient Ratio Tables (rji/rj'j): 

Having created the different tables for each antenna, we now create a 
new set of tables, which contain the ratios of the received signal amplitudes for 
the antenna pairs. For three receiving antennas we get two tables, while for N 
receiving antennas there will be — 1 tables. 

Figure 5.2 shows graphs of the ratio table data for sea state zero, 
at 9 GHz and a transmitter height of 20 meters. Note that for clarity only one 
“slice” is shown, referred to the particular transmitter height. The table is actually 
three-dimensional (range, transmitter height, and amplitude ratio). 

d. Amplitude Measurement: 

The first and major assumption is that all the receiving channels are 
calibrated. In this case, two similar signal amplitudes at the input to different 
channels will create the same amplitude at the outputs of those channels. 

Actually, there is no need for single channel amplitude measurements, 
since we deal with ratios. For this reason, we measure the individual signal ampli- 
tudes and immediately divide between pairs to form the desired ratios. 

e. Comparison: 

We now have in hand all the details needed to find the correct location 
of the transmitter, assuming of course, that the model is accurate. 

A comparison is made among the ratio levels from steps c and d. Every 
place that both levels are the same, is “marked down” by a “candidate location 
flag.” Figure 5.3 shows the received signal ratios (constant level line) and the 
corresponding candidate location flags, for four antenna pairs. 

Again, we should remember that the real comparison is being made 
between many more points, since there is a third dimension (height). For simplicity, 
we show here only one of 36 height slices. 
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RX. ANT. HEIGHTS; 14/12 METER RX. ANT. HEIGHTS: 16/12 METER 
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Figure 5.2: Attenuation Coefficient Ratio. 
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The dots in Figure 5.3 (upper portion) represent the real calculated 
values, while the connecting curves are made by extrapolation. Therefore, a can- 
didate location flag is set only when the measured signal ratio is within a specified 
distance from the calculated signal ratio. This distance is determined by the width 
of a comparison window which is a softvraxe parameter. The wider the window, the 
more candidate location flags we will set. On the other hand, we cannot allow too 
narrow a window, since there is a risk of loosing correct locations. It is obvious that 
as we use higher resolution in the lookup table, we can also narrow the comparison 
window. 

Figure 5.4 shows the candidate location flags for the same scenario as 
in Figure 5.3, but with a wider comparison window. Notice that more flags appear 
in this case. 

For a higher sea state, there are random fluctuations in signal level, 
and the comparison becomes harder. We must open the window wider, but at the 
same time more incorrect locations are flagged. A practical compromise is required, 
f. Validation: 

After we have all of the candidate location flags, we need to decide 
which location is correct. The method of doing this is by doing a cross location: 
comparing between all the flags, and checking which one repeats in all of the ratios. 

If after this process there is more than a single location (more than 
only one cross location), the measurement cannot give a correct location, so we 
discard it, and start a new session. 

Practically, using three receiving antennas (two ratios) gives no single 
location in most cases, and therefore we need to add more antennas in order to have 
more ratios. During the experiments we tried using different numbers of antennas, 
and the final sufficient result uses five receiving antennas. 
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Figure 5.5 shows the candidate location flags for each antenna pair and 
the crossing between them, giving a single valid location marked by a diamond. The 
sea state is zero, and the window width is 2.5%. Here, we can see not only a single 
“slice”, but the whole combination of range and height. From this it is obvious that 
a single pair of antennas gives no information at all (too many candidate location 
flags). 

For the case of three antennas, (see Figure 5.6), the ambiguity is re- 
duced significantly (only two flags overlapping), while adding a fourth antenna, 
reduces the location to a single, correct one. In this particular case, four receiving 
antennas give the correct location. In many cases, though, four antennas will still 
leave some order of ambiguity, but the fifth one will, in most cases, give a correct, 
single location. 

The ambiguity level reduction as a function of nximber of receiving 
antennas is well illustrated in Figure 5.6, for sea state zero. 

For a wider window there will be more flags, as shown in Figure 5.7, 
and here we can see that only the fifth antenna solves the ambiguity. 

When the sea state is higher, a wide window will lead to a large level 
of ambiguity. Figure 5.8 shows a 12% window in sea state 3, where there axe four 
locations. Reducing the width to 5% will reduce the number of locations to two, as 
shown in Figure 5.9, while the ambiguity will totally disappear for a 2.5% window, 
as in Figure 5.10. 

g. Output: 

The previous step brought up a single location flag, which shows the 
correct location of the transmitter in the range-height table. The output step 
translates the table location into the real coordinates of range and height, and 
types it as an output result. 
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C. FINAL RESULTS 



As mentioned earlier in this chapter, the final resolution used here, is 130 steps 
of 50 meters each, from the minimum range up to 7.5 Km, and then 125 steps of 
100 meters each to the maximum range. The comparison window width is 2.5%. 

With these parameters, estimation results are given in Figures 5.11 and 5.12, 
for sea states 0 and 3, respectively. It should be mentioned that the estimated 
range sometimes deviates from the actual range by a few tens of meters, which 
cemnot be seen in the graph, due to resolution. 
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RX. ANT. HEIGHT; 14/12 METER RX. ANT. HEIGHT; 16/12 METER 





Figure 5.3: Comparison Process with 2.5% Window. 
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Figure 5.3 (continued); Comparison Process with 2.5% Window. 
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RX. ANT. HEIGHT: 14/12 METER RX. ANT. HEIGHT: 16/12 METER 
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RX. ANT. HEIGHT: 18/12 METER RX. ANT. HEIGHT: 20/12 METER 






Figure 5.4 



(continued): Comparison Process with 12% Window. 
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CANDIDATE LOCATION FLAGS FOR 2.5 PERCENT WINDOW 

RX. ANT, HEIGHT: 14/12 METER RX. ANT. HEIGHT: 16/12 METER 
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Figure 5.5: Validation Process. Correct Location Marked by a 
Diamond. 
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CANDIDATE LOCATION FLAGS FOR 2.5 PERCENT WINDOW 

2 RECEIVING ANTENNAS 3 RECEIVING ANTENNAS 
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Figure 5.6: Ambiguity Level vs. Number of Receiving Antennas; Sea 
State 0 and 2.5% Window. 



57 



(KM) RANGE (KM) 



VALIDATED LOCATION FLAGS FOR 12 PERCENT WINDOW 
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Figure 5.7: Ambiguity Level vs. Number of Receiving Antennas; Sea 
State 0 and 12% Window. 
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Figure 5.8: Ambiguity Level vs. Number of Receiving Antennas; Sea 
State 3 and 12% Window. 
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CANDIDATE LOCATION FLAGS FOR 5 PERCENT WINDOW 

2 RECEIVING ANTENNAS 3 RECEIVING ANTENNAS 




OC 02 

(y3T3»N) iHOQH ‘Xy 



OC 02 01 

(yiON) iH3l3H ‘xy 




O 

CN 




u 

O 

z 

< 

a: 






ov oe oz OL 

IHOGH 'XH 



in 

< 



o 

2 

> 

LJ 

CJ 

UJ 

cr 



T 1 1 1 1 r 



J I I ! L 



o 

CN 



2 



UJ 

o 
2 
< 
GO cr 



OV OC QZ OL 

(d313n) iH9G!-l *Xd 



Figure 5,9: 



Ambiguity Level vs- Number of Receiving Antennas: 
State 3 and 5% Window. 
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VALIDATED LOCATION FLAGS FOR 2.5 PERCENT WINDOW 
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Figure 5.10: Ambiguity Level vs. Number of Receiving Antennas; Sea 
State 3 and 2.5% Window. 
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Figure 5.11: Estimated Results (Range); Sea State 0. 
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ESTIMATED VERSUS ACTUAL HEIGHT - SEA STATE 0 




Figure 5.11 (continued): Estimated Results (Height); Sea State 0. 
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Figure 5.12: Estimated Results (Range); Sea State 3. 
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Figure 5.12 (continued): Estimated Results (Height); Sea State 3. 
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VI. ADDITIONAL PARAMETERS CONSIDERATION 



This chapter deals with some parameters that have influence on the calcula- 
tions cind terms used for the estimation process. A few of them were mentioned 
but not discussed for all cases, and the others were not discussed previously. 

A. VERTICAL POLARIZATION 

So far we dealt with horizontal polarization, and all the calculations were done 
for this polarization. Vertical polarization complicates the problem a bit, since the 
reflection coefficient, p, changes significantly with the grazing angle, as shown in 
Figure 3.4. 

Let us check some representative cases in order to see how much the reflection 
coefficient really changes due to vertical polarization. 

Case 1 : transmitter is located 10 Km apart from the receiver, at a height of 
10 m, and it transmits with a frequency of 10 GHz. The receiver antennas are at 
heights of 12 m and 20 m. So, our case has the parameters; 
i2 = 10 km, 
hf = 10 m, 

/ = 10 GHz (A = 3 cm), 
hri = 12 m, and, 
hr 2 = 20 m. 

Using equations (5.4) to (5.12), we find that the grazing angles are equal to: 
V'l = 0.126°, and V’? = 0.172°, respectively for the receiving antenna heights. 
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Using the correct graph, (Figure 3.4), we find that the reflection coefficient due to 
the reflection, equals 0.98 for both grazing angles. 

Case 2 : in the second case the transmitter is higher, and the range between 
the receiver to the transmitter is reduced: 
i? = 4 km, 
ht = 40 m, 

/ = 10GHz (A = 3 cm), 
hri = 12 m, and, 
hr 2 = 20 m. 

Using the same method as in the first case, we find here, that the reflection 
coefficient for the lower receiving antenna is 0.82, while for the higher antenna it 
is 0.80. 

Case 3 : here we take an extreme case of very high transmitter, as well as 
receiver antennas: 

R = 4 km, 
ht = 100 m, 

/=10GHz (A = 3 cm), 
hri = 38 m, and, 
hr 2 = 50 m. 

In this case, the reflection coefficients are 0.58 for the 38 m receiving antenna, 
and 0.55 for the 50 m antenna. 

We can conclude that for low geometry (transmitter and receiver), and large 
ranges, the reflection coefficient cein be assumed the same for all receiving antennas. 
For shorter ranges and larger heights, we should take into account the changes in 
the reflection coefficient for the different receiving antennas. 
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Practically, for the vertical polarization case, we need to modify the off line 
tables by the amount of the reflection coefficient changes. This can be done, since 
we know how to calculate it with the given parameters, and for the whole range 
and transmitter height field. The only parameter that has doubt is the kind of 
polarization. This verification must be done by another system, and provided as 
an input to the estimation process. 

B. RECEIVER HEIGHT 

As we saw in Chapter 4, the higher the receiver antenna is, the more complex 
the curve is. In our work we used a maximum receiving antenna height of 20 meters, 
since it is a reasonable height for a conventional missile boat. For a larger ship (such 
as a corvette or destroyer, for example), we shall need higher receiving antennas. 
The only influence, as mentioned before, is on the curves complexity. The solution 
for this is to use a greater resolution in the look up tables, to compensate for the 
higher ambiguity. 

C. SHIP TILTING 

For the ships of our interest (relatively small missile boats), the worst case 
rolling angle can be around 15° for cross wind cruising direction in sea state 3. 
For sea state 4, the extreme angle is about 20° for the same conditions. Since 
the tilt occurs for both sides, we can say that, to a first degree of approximation, 
the average tilt is 0°. Even though, let us check the worst case. The geometry of 
the problem is given in Figure 6.1. We will now check the tilt impact on the look 
up tables for two cases: the closest antennas (12 meters and 14 meters), and the 
extreme antennas (12 meters and 20 meters). 
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Figure 6.1: Tilted Ship Geometry. 

Referring to Figure 6.1, in our case, hi = 12 m, /12 = 14 and /is = 20 m. 
The maximum tilt (roll) angle, On»ax> is taJcen as 15° for sea state 3, and 20° for 
sea state 4. 

From Figure 6.1 we see that the actual height of each receiving antenna, after 
being tilted, is: 

/i(- = /i,- • COSO • (i = 1, 2, • • • , 5). ( 6 - 1 ) 

For sea state 3 (Omax = 15°), we obtain; 

h'l = 12 • COS 15° — 11.6 m (where hi — 12 m), (6.2) 

and 

/15 = 20 • COS 15° = 19.3 m (where /12 = 20 m). (6-3) 

The actual distance between the antennas is 7.7 meters, instead of a theoretical 
value of 8 m. Figure 6.2 shows the ratio ^ 5 / 1/1 for receiving antennas located 
perpendicular to the sea level, in heights of 20 meter and 12 meter, respectively. 
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RX. ANT. HEIGHT: 20/12 METER 




RX. ANT. HEIGHT: 19.3/11.6 METER 




Figure 6.2: Ratio Between Attenuation Coefficients; Sea State 3, 
Receiving Antennas at 12 m and 20 m. 
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On the same figure we see also {tjI )' for the tilted case. Figure 6.3 shows the 
ratio between those ratios: 




Figure 6.3: Ratio Between Ratios of Attenuation Coefficients; Sea State 
3, Receiving Antennas at 12 m and 20 m. 

Now we will repeat the calculations for receiving antennas at heights of 14 
meters and 12 meters. Note that h[ remains the same (11.6 m), while h '2 becomes: 

h '2 = 14 ■ cos 15° = 13.5 m (where /i 2 = 14 m). (b-4) 
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Figure 6.4: Ratio Between Ratios of Attenuation Coefficients; Sea State 
3, Receiving Antennas at 12 m and 14 m. 

The actual distance here is 1.9 meters, instead of 2 meters. Figure 6.4 shows 
the ratio between r/l/r/i (perpendicular emtennas), and (vlyiiVi)' (tilted anten- 
nas). 

Repeating the same calculations for sea state 4 (omax = 20°), we obtain: 

h[ = 12 • cos 20° = 11.3 m (where hi = 12 m), (6.5) 

and 

/ig = 20 • cos 20° = 18.8 m (where hs = 20 m). (6-6) 
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The actual distance between the antennas is 7.5 meters, instead of a theoretical 
value of 8 meters. Figure 6.5 shows both ratios of t/I/t/i and while 

Figure 6.6 shows the ratio between those two ratios. 

For the closest antennas we obtain: 

/i 2 = 14 • cos 20“ = 13.2 m (where /12 = 14 m), (6.7) 

while h\ remains the same, 11.3 m. The actual difference is, therefore, 1.9 meters, 
instead of 2 meters. Figure 6.7 shows the ratio between (perpendicular 

antennas), and (^0^/(^i)^ (tilted antennas). 

To summarize these results, we see that the method is sensitive to ship tilting, 
since there is a difference between the tables for perpendicular and tilted config- 
urations. The conclusion of these results is that the receiving antennas must be 
stabilized. 

D. RECEIVER NOISE 

In the case of our work, the received signal to noise ratio was assumed to be 
relatively high. The ranges of interest axe less than 20 Km, while the missile ERP 
is in the order of 100 to 120 dBm. In addition, the transmitted signal propagates 
only one way. Therefore, any received signal is much larger than the receiver noise, 
and the signal to noise ratio is high enough, so we can disregard the noise. 

E. ATMOSPHERIC ATTENUATION 

The atmospheric attenuation can be treated in two ways: modify the look up 
tables with the attenuation data, since all the relevant parameters are known, or, 
more practically, to disregard it. Since we deal with ratios, the attenuation will 
have an influence only if the path difference is significant. As we already discussed 
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RX. ANT. HEIGHT: 20/12 METER 




RX. ANT. HEIGHT: 18.8/11.3 METER 




Figure 6.5: Ratio Between Attenuation Coefficients; Sea State 4, 
Receiving Antennas at 12 m and 20 m. 
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Figure 6.6: Ratio Between Ratios of Attenuation Coefficients; Sea State 
4, Receiving Antennas at 12 m and 20 m. 

in Chapter III, the path difference is negligible compared to the total path length. 
Therefore, the attenuation is practically the same for the direct and indirect path, 
so it has no influence on the problem, hence can be neglected. 

F. OTHER MODELS 

Appendix A gives other models for reflection coefficients, instead of the Ament 
model, that we used in this work. To use another model, the only thing that should 
be taken care of is modifying the model equation in the computer program. This 
point will be discussed in detail in the next chapter. 
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Figure 6.7: Ratio Between Attenuation Coefficients; Sea State 4, 
Receiving Antennas at 12 m and 14 m. 

G. EARTH CURVATURE 

In our work, we assumed a plane earth. In general, some geometries need 
to take account of the earth’s curvature. Our particular geometry results in no 
influence from the earth curvature, as shown below. The extreme condition is for 
the lowest receiving antenna, and the lowest missile. 

The maximum range that still allows plane earth assumption is calculated by 
using equation 6.8 [Ref. 7]. The problem geometry is shown in Figure 6.8. 

7?max = \J‘^0,ehr + + \J 2af^ht + , (b.8) 
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Figure 6.8: Plane Earth Assumption Geometry. 

where: 

Oe is the effective earth radius (= 4/3 x Re)^ 

Re is the earth radius (= 6370 Km), 
hr is the receiving antenna height, and, 
ht is the transmitter height. 

For all practical purposes, the term of under the square root can be ne- 
glected. 

Our extreme case is when the missile cruises in /i( = 5 m, and the lowest 
receiving antenna is located at height of ht = 12 m. For this geometry we obtain 
a maximum range of: 

i?max = \/2 • 8493 • 10^ • 12 + \/2 • 8493 • 10^ • 5 = 23.5 Km. 
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Since we deal with relative ranges of up to 20 Km, we can disregard the earth 



curvature. 
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VII. COMPUTER PROGRAM 



The computer program in this work was written in APL. The listing is given in 
Appendix B and it includes all the necessary comments. This chapter presents an 
overview of the whole program, and then explains each function in greater detail. 

The program consists of a few functions which call each other in sequence, and 
a main management function. 

The sequence follows exactly after the one described in Chapter V, under 
“Detailed Discussion.” All the inputs that must be given by external means, are 
inserted here by the user, as software parameters. 

The program is started by typing RUN. The user is asked to enter sea state 
(0 to 7), and with this data, the RUN function calls the main function for receiving 
data, named DATAIN. This function automatically creates more needed parame- 
ters, in a random manner, between given limits. It is also possible for the user to 
dictate each one of these parameters, instead of creating them randomly by the 
program: 

• frequency - between 2 to 18 GHz, 

• ERP (Effective Radiated Power) - between 100 to 120 dBm (10^ to 10^ 
Watt), 

• initial range for simulation (range where missile starts to radiate) - between 
10 to 20 Kilometers, and, 

• transmitter height - between 5 to 40 meters. 

Other relevant parameters are programmed off-line as software parameters, 
and can be changed easily at any time. Those are: 
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• window width, as discussed in Chapter V, 

• minimum trajismitter height for the look up table, 

• number of steps for the transmitter height for the table, 

• resolution of transmitter height steps for the table, 

• minimum range for simulation (if the missile passed this range without being 
intercepted, the simulation is failed and stops), 

• range resolution for the look up table, both high resolution for short range, 
and regular resolution for longer range, 

• number of range steps for each range resolution, 

• resolution of simulation steps (range distance between each receiving), and, 

• locations (heights) of receiving antennas. 

The function then calculates the maocimum range, and shows the major pa- 
rameters on the screen (a full session is demonstrated later in this chapter.) 

With all the necessary data, DATAIN returns to the management function, 
RUN. In the next step, RUN sends all the relevant parameters to the main function 
STARTPRO, and activates it. 

STARTPRO is responsible for all the major calculations: 

• creating off line tables which include also the attenuation coefficient, rj^ , 
and the reflection coefficient, p, 

• receiving five signals in the receiving antennas, and divide their powers to 
determine the ratios, 

• comparing between the off line tables and the actual ratios, 

• marking the correct (single) location, and finally, 

• sending the location to the management function. 
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The signals for STARTPRO are created by the simulator SIM. This function 
is activated by STARTPRO, every time it needs a new set of receiving signals. 

SIM receives the relevant parameters from DATAIN, and creates five signals 
due to those parameters and the geometry. The signal is calculated using equa- 
tion 7.1. 



5.- 



= ERP • r/2 . p . 



1 

47tR2 ’ 



(7.1) 



where: 

ERP is the effective radiated power of the transmitter, 

R is the relative range between the transmitter and the receiver, 

T}^ is the attenuation coefficient, calculated by using Equation (3.5), and, 

p is the reflection coefficient, calculated by using Equation (3.20). 

After the management function RUN received the location of the transmitter 
from STARTPRO, it calls the translating function RHFIND in order to translate 
the coordinates of the table into range and height, and to type it. 

A conceptual flow chart is given in Figure 7.1, and a simulation session is 
shown in Figure 7.2. 
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STARTPRO 



SIM 



RHFIND 



_\l^ 



Transmitter 
Range and 
Height 




Figure 7.1: Estimation Process; Conceptual Flow Chart. 
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*RUN 

ENTER SEA STATE: 

0 - CALM 1 - SMOOTH 2 - SLIGHT 

3 - MODERATE 4 - ROUGH 5 - VERY ROUGH 

6 - HIGH 7 - VERY HIGH 



*2 (for example) 

RANGE INTERVAL FOR SIMULATION: 

MINIMUM RANGE (IN METERS): 1000 

MAXIMUM RANGE(IN METERS): 20000 

SEA STATE: 

2 

FREQUENCY (IN GHZ): 

8.85 

INITIAL RANGE (IN KILOMETERS): 

14.5 

TRANSMITTER HEIGHT (IN METERS) - ACTUAL: 
36 

RECEIVER ANTENNAS HEIGHTS (IN METERS): 

12 14 16 18 20 

TRANSMITTER RANGE (IN KILOMETERS): 

14.4 

TRANSMITTER HEIGHT (IN METERS): 

36 

UPDATED ACTUAL RANGE (IN KILOMETERS): 
14.4 



Figure 7.2: Simulation Session Example. (Lines Starting with Asterisk 

(*) are Entered by the User.) 
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VIII. SUMMARY AND CONCLUSIONS 



The work which had been done during this thesis research, together with the 
results obtained, show that the suggested method for passive range estimation 
using multipath effect can probably be used. The work developed the method 
theoretically, and proved its feasibility. 

The accuracy of the estimated range and height depends upon many param- 
eters, as well as receiving antennas configuration. A trade off should be done be- 
tween number of receiving antennas, their heights, separation and window width, 
as described in Chapter V. 

There is no analytical way for finding the optimum trade off, and only trial 
and error sessions can achieve a local optimum, as was done during this research. 

Another question of interest is the accuracy of the rough sea model (Ament 
model in our case). 

The answer to this question can only be obtained by sea trials. Those trials 
should answer the questions listed below: 

• Is the theoretical analysis of the additional parameters in Chapter VI cor- 
rect, or should we modify these parameters? 

• What is the estimation accuracy if an incorrect sea state is entered, and 
is there a preferred direction for entering sea state (is it better to enter a 
higher sea state in case of uncertainty, or lower state)? 

• How does the wind direction influence the accuracy, if at all, (upwind, down- 
wind, and crosswind)? 

• What is the practical feasibility of the suggested method, and with which 
limitations (such as sea state and wind direction)? 
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We should notice that the Ament model neglects the shadowing effect, al- 
though it deals with small grazing angles, where this effect does occur. It is ques- 
tionable how far the shadowing effect changes the reality from the model. 

A possible influence of the shadowing effect is the fact that the amplitude of 
the reflected wave will be lower, so that the signal minima will not be so sharp as 
is shown in Figure 4.2. A decrease in the amplitude of the specular reflection will 
also occur occur because of high sea state, and capillary waves. If sea trials show 
consistent errors, we can conclude that the Ament model is inadequate. In this 
case it will be necessary to develop another more accurate model. 

It can be predicted that for sea state 0, the results will be more accurate than 
in any other case, while the more roughness, the less accurate the results. 

In order to improve the estimation process, and to reach a lower rate of errors, 
we can modify the program. The modified program will consider also the past 
readings, and compare each present estimate with the previous ones. If the range 
becomes shorter by a reasonable increment, and the transmitter height remains 
almost the same, we can conclude that the estimate is correct. On the other hand, 
if the estimate does not meet this test, then another estimate should be made, and 
compared, and so on. 
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APPENDIX A: REFLECTION COEFFICIENT MODELS 



In Chapter III we discussed the most popular model for reflection coefficient 
of a rough sea, the Ament model. There are several other models, which are given 
below. It should be noted that there axe additional models which we will not cover 
here. 



A. BULLINGTON MODEL 



This model distinguishes between different surface profiles [Ref. 3,12]; 

1. Sawtooth Profile 

The wave height is distributed uniformly, but does not necessarily have 
uniform slope. The phase values have equal probability between —(f>m and +<pmi 
where 



4^m 



2Tr<Th sintp 
A 



(A.l) 



ah is the RMS deviation of the waves height, 
is the grazing angle, and, 

A is the signal wavelength. 



The reflection coefficient for the sawtooth profile is given by; 



sin . 

P = —7 = sinc^m, 

0m 



and is shown in Figure A.l. 



(^. 2 ) 
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Figure A.l: Bullington Model — Sawtooth Profile. 

2. Sinusoidal Profile 

The reflection coefficient in this case is given by Equation (A. 3), and is 
shown in Figure A. 2. 

(A3) 

where Jq is the Bessel function of the first kind and zero order. 

3. Triangular Profile; 

Here, the reflection coefficient is given by: 



P = 



' sin<i)t 



(A.4) 
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Figure A. 2: Bullington Model — Sinusoidal Profile. 

where the maximum phase is: 



<i>t = 



7T<7fc sin xp 
A 



(A.5) 



Figures A. 3 and A. 4 show the triangular profile refiection coefficient. 

4. Sine Square Profile 

The refiection coefficient is more complex in this model, and is given by 
Equation (A. 6), and shown in Figure A. 5. 



p = cos 4>t ■ Jo{(f>t)j 



(A.6) 



where <j>t is the same as given in Equation (A. 5). 
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Figure A. 3: Bullington Model — Triangular Profile. 



B. DAVIES MODEL 



The Davies model assumes that the surface has a random behavior, which 
yields a gaussian height distribution and correlation coefficient. This model neglects 
the shadowing effect, and is good for small to medium grazing angles (0 < 40°). In 
order to describe the surface, we introduce an arbitrary parameter, m/a. Experi- 
ments above the sea surface show that a value of m/a = 50 to 60, gives sufficient 
correlation between theoreticcil and experimental results. [Ref. 3]. 

The reflection coefficient for this model is given by: 



P = 



m 

2(7 



y/sectp 



exp 




(.4.7) 



Figure A. 6 shows the reflection coefficient. The value m/2a is 27 {m/a = 54). 
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Figure A.4: Bullington Model — Triangular Profile (expanded). 

C. ISAKOVICH MODEL 

In this model we assume that the surface is irregulcir, isotropic in all directions, 

and heis a normal height distribution. The parameter m/<7 describes the surface (as 

in the Davies model). It is the only parameter which dictates the intensity of the 

scattering in a given direction. Changing the ratio will change the profile shape, 

as shown in Figure A. 7 (m/cr = 10), and Figure A. 8 (m/cr = 20). The shadowing 

effect is neglected and only far-field is calculated. The Isakovich model calculates 

the dispersion of waves for random surface [Ref. 3]: 

1 — cos V’ cos Et cos 9 + sin t/> sin Et 
(sin V* + sin Ft)^ 

1 cos^ Ip + cos* Et — 2 cos ip cos Et cos 9 

^ _8 V <7 / (sin 0 + sin Ff )* 



(A.8) 
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Figure A. 5: Bullington Model — Sine Square Profile. 

This model is much more complicated than the others, since it has two more 
parameters: 

Et is the observation angle, and, 

9 is the azimuth angle. 



In order to illustrate results comparable to those for the other modes, we 
assume that the azimuth angle is zero {9 = 0). For simplicity, we calculate the 
reflection coefficient for two cases: zero observation angle, Et = 0° (transmitter 
and receiver at the same level), and Et = —30® (transmitter higher than the 
receiver). 
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Figure A.6: Davies Model. 

Using the assumption ^ = 0, we simplify Equation (A. 8) to; 
1 — cos xj) cos Et + sin ip sin Et 



m 

P»=o> = — 



(sinip + sin Et^ 



exp 



fi fniV. 


cos^xp + cos^ Et — 2 cos xp cos Et 


[s V(7 J 


(sin 0 + sinE't)^ 



(. 4 . 9 ) 



For the case of zero observation angle {Et = 0°), Equation (A. 9) becomes even 
more simple: 

m / 1 — cos ip\ cos^ 0 + 1 — 2 cos xp 



P E. =0« 



<7 V sin^ xp 






sin^ xp 



(. 4 . 10 ) 
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This case is illustrated in Figures A. 7 and A. 8 for two different values of 
mfo. For —SO” observation angle Equation (A. 11) applies. The result is shown in 
Figure A. 9. 



P Et=-30‘> 



m 2 — y/3 cos rp — sin tp 
2(sint/’ — 0.5)^ 



• exp — 



1 /m\2 


cos^rl) + 0.75 — a/S cos V’ 


8 \~a) 


(sin t/» — 0.5)^ 



(All) 




Figure A. 7: Isakovich Model with Zero Observation Angle (Et = 0°), 
m/cr = 10. 
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REFLECTION COEFFICIENT 




GRAZING ANGLE (MILRAD) 



Figure A. 8: Isakovich Model with Zero Observation Angle, {Et = 0°), 
m/(7 = 20. 
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REFLECTION COEFFICIENT 



o 




GRAZING ANGLE (MILRAD) 



Figure A. 9: Isakovich Model with Observation Angle of — 30°,m/<r = 20. 
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APPENDIX B: COMPUTER PROGRAM LISTING 



This appendix gives the listing of the computer program described in Chapter 
VII. RUN is the management function, which calls the other functions in order. 
DATAIN is the function responsible for the input data for further calculations. 
STARTPRO is the main function which does all the major calculations, as well as 
creating the look-up tables, and comparing the ratios of the actual received signals 
with the calculated ratios. SIM is the simulation function, which simulates the 
received signals from the missile. RHFIND is the function which translates the 
location of the missile in the table into a range and height. 



**RUN 



CO] 


RUN 




Cl] 


0 MANAGMENT FUNCTION 




C2] 


'ENTER SEA STATE'.' 




C3] 


t 1 




C4] 


'Q-CALM 1-SMOOTH 2-SLIGHT' 


C5] 


^2 -MODERATE n -ROUGH 5 -VERT 


ROUGH 


C6] 


^6 -HIGH 7 -VERY HIGH^ 




C7] 


1 t 




C8] 


SEA^^U 




C9] 


p ENTERING ROUGHNESS PARAMETER 




CIO] 


t 1 




Cll] 


DATAIN 




C12] 


p INITIALIZATION 




C13] 


RR:F STARTPRO HRS 




C14] 


p STARTING 




->C15] 


^iFLAG=l)/STO 




C16] 


p END OF SESSION IF RANGE CANNOT BE 


FOUND 


C17] 


RHFIND 




CIS] 


p GO TO CHECK REAL RANGE AND HEIGHT 




-^Cl9] 


STOx^O 





96 



**DATAIi:^ 

[0] DAT AIN 

[1] R INITIALIZATION FUNCTION 

[2] fl □DODODDO 

[3] MINHEICT-^n 

[4] fl 'MINHEICT' IS THE MINIMUM TX, HEIGHT (LESS 1 METER) 

[5] fl □□□□□□□□ 

[6] HEIGTST-^ZS 

[7] 0 'HEICTST' IS THE NUMBER OF HEIGHT STEPS 

C8] fl □□□□OODO 

[9] HEIGTRES-^1 

[10] fl 'HEIGTRES' IS THE RESOLUTION OF HEIGHT STEPS 

[11] fl □□□□□□□□ 

[12] C-«-0.975 

[13] 0 □□□□□□□□ 

[14] Dt-1.025 

[15] fl ARE LOWER AND UPPER LIMITS FOR OK (STARTPRO) 

[16] fl □□□□□□□□ 

[17] RMIN^IOQO 

[18] n 'RMIN' IS MINIMUM RANGE 

[19] RMINS^RMIN 

[20] fl 'RMIN S' IS MINIMUM RANGE FOR SIMULATION 

[21] fl □□□□□□□□ 

[22] NSTEPSf-125 

[23] fl 'N STEPS' IS THE NUMBER OF RANGE STEPS 

[24] fl □□□□□□□□ 

[25] RESO-<-100 

[26] fl 'RESO' IS THE RANGE STEPS RESOLUTION (LENGTH OF EACH STEP) 

[27] fl □□□□□□□□ 

[28] HNSTEPS-<rlZ0 

[29] fl 'HNSTEPS' IS THE NUMBER OP HIGH RESOLUTION RANGE STEPS 

[30] fl □□□□□□□□ 

[31] HRESO^SQ 

[32] fl 'HRESO' IS THE RANGE STEPS (HIGH) RESOLUTION FOR SMALL 

[33] n RANGES 

[34] fl □□□□□□□□ 

[3 5] SIMSTEP->-50 

[36] n 'SIMSTEP' IS THE RANGE STEPS RESOLUTION IN THE SIMULATION 

[37] n FUNCTION (VSIM) 

[38] F-e-(1999+(?17001))Tl000 

[39] fl □□□□□□□□ 

[40] flF-«-9.35 

[41] fl <F' IS FREQUENCY (IN GHZ)’, 2rl8 GHZ 

[42] EFP-e-lOOOOOOOxPlOO 

[43] fl □□□□□□□□ 

[44] flEi?P-e-100000000 

[4 5] fl 'ERP' IS TRANSMITTER EFFECTIVE RADIATED POWER (IN WATTS)', 

[46] fl 1E741P9 WATT (100^120 DBM) 

[47] PJ-*-20000x(500+l?300)^1000 

[48] fl 'RI' IS INITIAL RANGE (IN METERS); 10,000*16,000 METER 

[49] 0 □□□□□□□□ 
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□□□□□□□□ 



C50] fi/?I-*-15400 

[51] HAT*-^ + ?36 

[52] 0 

[53] oHAT-<r22 

[54] Q 'HAT' IS TRANSMITTER HEIGHT (.IN METERS); 5t4Q METER 

[55] o/fi?3-*-14 + 5?16 

[56] 0 'HR3' ARE RECEIVER'S ANTENNAS HEIGHTS (IN METERS) 

[57] 0 'HR3' IS BETWEEN 15^30 METERS 

[58] fi □□□□□□□□ 

[59] HR3->- 12 14 16 18 20 

[60] 'RANGE INTERVAL FOR SIMULATION;' 

[61] 'MINIMUM RANGE (IN METERS): ' ,^RMIN 

[62] 'MAXIMUM RANGE (IN METERS): 

' , !sRMIN-i- (NSTEPS>(RESO ) + (HNSTEPS^HRESO ) 

[63] 'SEA STATE:' 

[64] SEA 

[65] fi 'SEA' IS ROUGHNESS PARAMETER 

[66] 'FREQUENCY (IN GHZ):' 

[67] F 

[68] 'INITIAL RANGE (IN KILOMETERS):' 

[69] i?I-1000 

[70] 'TRANSMITTER HEIGHT (IN METERS)-ACTUAL:' 

[71] HAT 

[72] 'RECEIVER ANTENNAS HEIGHTS (IN METERS):' 

[73] HR3imR3l 

[74] p DISPLAING THE INITIAL SIMULATION DATA 

■*■[75] -*-0 
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**STARTPRO 

CO] FO STARTPRO HR;DIO ;HA;R; LAMBDA ;ET1 ;ET2 ; ET3 ; ETn -,ET5 ;ET12 ‘.ETlS ; 

Fri4;Fri5;SI;SJ12;SI13;SI14;SI15;Fl2;Fl3 :Fl4;i?15;OA'12; 
OK13;OKin; OKI 5 ; TEMP ; TEMPI ; TEMP2 ; SH RDl ; RD2 ; RD3 ; FD4 ; 

RD 5 ; DELTA 1 ; FFLM 2 ; FEiM 3 ; OFEM 4 ; DELM 5 ; TETA 1 ; TETA 2 ; 
TETA3iTETAK; TETA 5 ; XII ;XI2 ;XI3 ; XI^ ’,XI5 xRROl •,RR02 ; RH03 ; 
RH0'i;RH05 

Cl] fl MAIN CALCULATION AND COMPUTATION FUNCTION 

C2] QIO^l 

C3] R-^RMIN+(HRESOx\HNSTEPS) , (.HRESO^HNSTEPS )+RESOx xNSTEPS 

C4] n 'R< IS RANGE VECTOR ('HN STEPS' STEPS OF 'HRESO' METER 

C5] n EylCE, PLUS 'N STEPS' STEPS OF 'RESO' METER EACH), 

C6] n STARTING AT 'RMIN' 

C 7 ] HA-f-MINHEIGT+HEIGTRES^ x HEIGTST 

C8] fi 'HA' IS TRANSMITTER HEIGHT VECTOR ('HEIGTST' STEPS OF 
C9] n 'HEIGTRES' METER EACH) 

CIO] GRH*-(HEIGTST , (NSTEPS-^HNSTEPS ) ) p 1 

Cll] n 'GRH' IS MATRIX WHICH SHOWS ALL THE LOCATIONS OF 

Cl 2] R CANDIDATE RANGE AND HEIGHT COMBINATIONS 

C13] LAMBDA*-Q.2iF0 

C14] fl 

CIS] RDl-^iUHA-HRll) )*2)o.+F*2)*0.S 

C16] RD2*-(((HA-HRl2) )*2 )« . +A*2 )*0 . 5 

C17] RD3-fr(((HA-HRi3') )*2 ) « . +F*2 )*0 . 5 

C18] ED4-f-(((EA-EEC4] )*2 )« . +A*2 )*0 . 5 

C19] RD3ir(((EA-HRL3) )*2 )o . +E*2 )*0 . 5 

C20] fl 'RD' IS THE DIRECT DISTANCE BETWEEN TRANSMITTER 

C21] fl AND RECEIVING ANTENNAS (RESPECTIVELY) 

C22] Z)EE2’ill-*-(2xMxAECl] )® .tE 

C23] DELTA2-^(2^HAy-HRl2l)o .iR 

C24] DEE2’43-«-(2xEAxEEC 3] )o.+E 

C25] DELTA'i*-(2y-HAy-HRL'M)o .iR 

C26] DELTA3*(2y.HAy.HRZS))^ .iR 

C27] fl 'DELTA' IS THE DIFFERENCE BETWEEN THE DIRECT 

C28] fl PATH AND THE REFLECTED PATH FOR EACH RX. 

C29] fl ANTENNA (RESPECTIVELY) 

C30] Z’ET41-«-"lo(iS!(<}>pEDl)p |EECl]-EiO+ECl 

C31] Z’Er>l2-t-"lo(iS! (<t>pEC2 )p | EEC2] -HA)iRD2 

C3 2] rETA3-t-"lo((9(<})pED3)p|EEC3]-EiO*ED3 

C33] Z’EZ’A4-*-”lo(l5(<|)pED4)p |EEC4]-EiOTED4 

C3 4] rEri55-e-"lo(IS!(<})pED5)p|EEC5]-H.4)*EE5 

C35] fl 'TETA' IS THE OBSERVATION ANGLE BETWEEN THE RX. 

C36] fl ANTENNA AND THE TRANSMITTER (FROM THE LATEST 

C37] fl POINT OF VIEW)- FOR EACH RX. ANT. (RESPECTIVELY) 

C38] XJl-f-(i5(<j>pEDl )pEECl] +HA )*RD1+DELTA1 

C39] XI2-^(>st(<^pRD2)pHRl2']-i-HA)*RD2*DELTA2 

C40] Xl3-*-(i5(<t)pED3 )pHRL3)+HA)rRD3+DELTA3 

C41] Xl4-c(i?(4)pEE4 )pHRW+HA)^RD'i+DELTAn 

C42] XI5-e-(S)((t>pEE5 )pHRL51+HA )iRD5+DELTA5 

C43] fl 'XI' IS THE SINE OF THE GRAZING ANGLE FOR EACH 

C44] fl RECEIVING ANTENNA (RESPECTIVELY) 
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^[45] 


^(.SEA^O )/S0 








^[46] 


-(SEA=1 )/Sl 








^[47] 


->(SEi5 = 2)/S2 








^[48] 


->(SEA = 3 )/S3 








^[49] 


->(SE4 = 4)/S4 








^C50] 


->(S£’i9 = 5 )/S5 








^C51] 


->(SEi? = 6 )/S6 








^C52] 


*(.SEA = 7 )/S7 








C53] 


P 








[54] 


S0:SH^SHL^SHH^0 








^[55] 


■^CONT 








[56] 


Sl:Sff->0.03 2 5 








[57] 


SHL^O 








[58] 


SHH^O.065 








4-[59] 


■*CONT 








[60] 


127S 








[61] 


SHLf-O.O&S 








[62] 


21 








^[63] 


■*CONT 








[64] 


S3 :SH-«-0.26 5 








[65] 


SHL^O.21 








[66] 


SHH<’0.32 








^[67] 


^CONT 








[68] 


S4:SH^0.43 








[69] 


SHL^O.32 








[70] 


SHH^O.5^ 








->[71] 


^CONT 








[72] 


S5:SH^0.7 








[73] 


SHL^Q.S^ 








[74] 


SHH^O.Se 








->[75] 


^CONT 








[76] 


SQiSH^l.OB 








[77] 


SBL^O.Se 








[78] 


SHH^l.3 








^[79] 


’^CONT 








[80] 


57:5^^1.75 








[81] 


SHL^l . 3 








[82] 


SHH^2 . 2 








->[83-] 


^CONT 








[84] 


p ^SEL^ AND 'SEH^ ARE THE 


LOWEST 


AND 


' HIGHEST 


[85] 


fi (.RESPECTIVELY) RMS WAVES HEIGHT FOR 


EACH 


SEA STATE 


[86] 


CONT: 








[87] 


RH01^*~2^((o2)^SHxXI1*LAMBDA)*2 








[88] 


RH02^*~2^ ( (o2 )>(SH-xXI2iLAMBDA)*2 








[89] 


RH03-(-*~2x((o2)xSHxXI3*LAMBDA)*2 








[90] 


Hff04-*-*~2x( (o2 )^SH>^XiniLAMBDA )*2 








[91] 


RH03-<r*~2y.((o2)y-SH>.XISiLAMBDA)*2 








[92] 


p 'RHO' IS THE REFLECTION COEFFICIENT 


FOR 


EACH 


[93] 


p RX. ANTENNA (RESPECTIVELY) - 


CALCULATED AFTER 


[94] 


p AMENT MODEL (SEE VSIM) 
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[95] 

C96] 

C97] 

C98] 

C99] 

ClOO] 

ClOl] 

C102] 

C103] 

[104] 

[105] 

[106] 

[107] 

[108] 

[109] 

[ 110 ] 
[ 111 ] 
[ 112 ] 

[113] 

[114] 

[115] 

[116] 

[117] 

[118] 

[119] 

[ 120 ] 
[ 121 ] 
[ 122 ] 

[123] 

[124] 

[125] 

[126] 

[127] 

[128] 

[129] 

[130] 

[131] 

[132] 

[133] 

[134] 

[135] 

[136] 

[137] 
>■[138] 

[139] 

[140] 
"[141] 
"[142] 

[143] 

[144] 

[145] 

[146] 
"[147] 



rUAMBDA^R)) )*2 
r(.LAMBDAxR)))*2 
t(.LAMBDA^R)))*2 
t{LAMBDA^R) ) )*2 
tUAMBDA^R)) )*2 



BETWEEN THE MULTIPATH 
ANTENNAS, RESPECTIVELI 



ETl-^RHOly(lo(((o2 )x(rHA*2)+HRZll^HA ) )o 
ET2^RH02x{loU(o2 )x({HA*2)+HRL21xHA ) )« 

5r3^/?W03x (io( ( (o2 )x( (Hyl*2 )+n’/?[3]xff^ ) )» 

£T4-;?i704x (io( ( (o2 )x((HA*2)+HRZk2xHA ) )<> 

ET5^RH05xilo(. ((o2 )x((M*2 )+H/?[5]xM ) )o 
p <ETl'-r'ET5' ARE MULTIPATH ATTENUATION FACTORS FOR 
p EACH RECEIVING ANTENNA, RESPECTIVELY (MATRIX) 
ET12*-ETlrET2 + lE~13 
ET13-<-ETli-ET3 + lE~13 
ETin-(-ETl-i-ET'i + lE~13 
ETl5-<-ETli-ET5 + lE~13 
p '£T12 ' T '£ri5 ' ARE RELATIONS 
p FACTORS FOR EACH PAIR OF RX. 

R (MATRIX ) 

RCVxSI*SIM 

R 'SI' ARE RECEIVED SIGNAL LEVELS FOR ALL THE RECEIVING 
p ANTENNAS (VECTOR) 

REL i 

SJ12-*-(pCra)pSJ[l]TSI[2]+l£:”13 
SJ13 -*-(pC;?H)pSJ[1]tSI[3] +1£”13 
SJ14-(-(pG;?H)pSI[1]tSI[4]+1E"13 
SJ154-(pC/?H)pSJ[1]tSI[5]+1£:"13 
R 'SJ12 ' ^ 'SI15 ' ARE RELATIONS BETWEEN RECEIVED SIGNALS 
R PAIRS (MATRICES ) 

COMP: 

R12-^SI12rETl2+lE~13 
R13-^SI13 tET13-HE~13 
H14-(-SJ14t£T 14+1£'"13 
/?15-«-SI15 + Eri5 + l£’”l3 
p ';?12'r'i?15' ARE RATIO MATRICES 

p AND THEORETICAL VALUES 

0K\2->r(Rl2'^C)f^(R12<.D) 

OKl3-<-(Rl3^C)f^(R13^D) 

0/i:14-(-(^14 2C)a(H142D) 

0/fl5f-(ffl5SC)A(H15£D) 

p '0X12 ' + 'OKI 5' ARE LOCATIONS OP RECEIVED 
p WHICH ARE EQUAL 

p (ACTUAL RELATIONS 

p (MATRICES) 

MAT COMP: 

Oi?Xf-OXl 2 aOXI 3 AOXl 4 AOXl 5 

P 'CRH' ARE LOCATIONS WHICH ARE EQUAL ALSO BETWEEN 
P ANTENNAS PAIRS 

■^(l--^/ , CRH)/ RHFIN 

p CHECK IF THERE IS A SINGLE COMBINATION OF RANGE 
R AND HEIGHT 

*(FLAC-1 )/RHFIN 
*RCV 



BETWEEN ACTUAL 



SIGNALS 

TO THE PRE-CALCULATED VALUES 
EQUAL TO THE THEORETICAL ONES) 



IF MORE THAN ONE COMBINATION , GO TO ANOTHER RECEIVING 
SESSION (AS LONG AS THE RANGE IS HIGHER THAN THE 
MINIMUM RANGE). IN CASE OF A SINGLE COMBINATION GO TO 
TRANSLATE INTO REAL READING (RANGE AND HEIGHT) 



P 
P 
P 
P 

RHFINi^Q 
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**SIM 

CO ] A*-Siyi ; LAMBDA ; A1 ■, A2 \ A2 \An • A5 \ SXI ; SHR : SRHO 

Cl] p RECEIVED SIGNAL AMPLITUDES SIMULATION UT EACH 

C2] p RECEIVING ANTENNA), INCLUDING TRANSMITTER 

C3] p MOVEMENT INFLUENCE 

C4] LAMBDA-^-0 .2iF 

C5] FLAG*-0 

C6] SXI^~2o{(.\HR2-HAT)iRI) 

C7] p 'SXI' IS THE GRAZING ANGLE 

->-Ca] -*(.SEA = 0 )/LAB 1 

C9] SHR-^SHL+(.l? (.SHH-SHL)xl000 )^1000 

CIO] p 'SHR' IS THE RMS NAVES HEIGHT FOR SIMULATION 

-►Cll] ^LAB2 

C12] LABl:SHR-^0 

C13] LAB2:SRHO^*~2y.((,(,o2 )xSHR^SXI t LAMBDA) )*2 

C14] P 

CIS] P THE SIMULATION USES AMMI MODEL FOR CALCULATING 
C16] P THE REFLECTION COEFFICIENT (p): 

C17] p 
CIS] p 

C19] P p =£XP. (-2x(2xPJxH xSIN , (.XI)^LAMDA)*2) 

C20] p S RMS 

C21] P 
C22] p 
C23] 

-a.l-«-SPflOCl] X (£PPx4x (io( ( (o2 )x (PP3 Cl] +HAT)xRAT)i- (.LAMEDAxRI ) ) )*2 )^- ( 
04 )xRI*2 

C24] 

A2-»-SPH0C2]x(PPPx4x(1o(((o2 )x (^RR3 Z2l +HAT )xHAT ) r (LAMBDAxRI ) ) )*2)i ( 
04 )xRI*2 

C25] 

A34-SPP0C3]x(PPPx4x(1o( ((02 )x(ffP3 Z3l+HAT)xHAT)^ (LAMBDAxRI)) )*2 )^( 
04 )xpj*2 

C26] 

An*-SRHOW X (PPPX4X (lo ( ( (o2 )x (PP3 C4] ■^HAT )xHAT) ^ (.LAMBDAxRI ) ) )*2 )^ ( 
04 )xpx*2 

C27] 

715-^-SPflOC5]x(PPPx4x(io(( (02 )x(HR3L51+HAT)xHAT)t(LAMBDAxRI)) )*2 )t( 
04)xPJ*2 

C28] p 'SRHO' IS THE 'REFLECTION COEFFICIENT' (p) FOR 

C29] p EACH RECEIVING ANTENNA (RESPECTIVELY) 

C30] /lfAl,il2,/13,/14,A5 

C31] p 'A' IS A VECTOR OF RECEIVED AMPLITUDES (IN EACH 

C32] P ANTENNA-RESPECTIVELY) 

C3 3] RI-<^RI-SIMSTEP 

C34] p UPDATING THE TRANSMITTER RANGE (IN ' SIMSTEP' METER STEPS) 
-C35] ■*(RI^RMINS)/RET 
C36] ' ' 

C37] 'END OF SIMULATION SESSION'.' 

C38] FLAG*-1 

C39] p RANGE HAD NOT BEEN FOUND YET, BUT IS LESS THAN THE 

C40] p MINIMUM RANGE, OR END OF SIMULATION SESSION 

-C41] RETi^O 



102 



**RHFIND 

C 0 ] REF IN D ; IN DU ; IN DR ; MATER ; IND VEC ; DOC ; V EC ORE 

Cl] R FINDING AND DISPLAYING TEE RANGE AND EE I GET 

[2] fl OF TEE TRANSMITTER 

C 3 ] INDE-<r 10 00 0 O^MINEEIGT+EEIGTRESx i EEIGTST 

[4] fl 'INDE' IS AN EE I GET INDEX 

C 5 ] INDR^RMIN + ( ERE SO x x EN STEPS ) , ( ERE SO x EN STEPS ) + (RE SO x i N STEPS ) 

[6] fl 'INDR' IS A RANGE INDEX 

C 7 ] MATER^INDE « . +INDR 

C8] fl 'MATER' IS A RANGE / EEIGET INDEX MATRIX 

[9] INDVEC-^ , MATER 

[10] FIND: 

[11] VECCRE*,GRE 

[12] fl CEANGING MATRICES INTO VECTORS 

[13] LOC->-VECCRE/INDVEC 

[14] fl SELECTING TEE CORRECT LOCATION OF TEE TRANSMITTER 

[15] fl ('LOC' IS TEE TX. LOCATION) 

[16] ;?-<-iooooo Uoc 

[17] S-e-LLOCT 10000 0 

[18] fl TRANSLATING INTO REAL RANGE ('R') AND EEIGET ('E') 

[19] ' ' 

[ 20 ] ' ' 

[21] 'TRANSMITTER RANGE (IN KILOMETERS):' 

[22] /?4-1000 

[23] 'TRANSMITTER EEIGET (IN METERS):' 

[24] E 

[25] " 

[26] ' ' 

[27] 'UPDATED ACTUAL RANGE (IN KILOMETERS):' 

[28] (SJMSrEPxlOOO )+i?J+1000 

-*•[29] -*-0 
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